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Abstract 

We present a method for constructing new families of solvable one-dimensional 
diffusions with linear drift and nonlinear diffusion coefficient functions, whose tran- 
sition densities are obtainable in analytically closed-form. Our approach is based on 
the so-called diffusion canonical transformation method that allows us to uncover new 
multiparameter diffusions that are mapped onto various simpler underlying diffusions. 
We give a simple rigorous boundary classification and characterization of the newly 
constructed processes with respect to the martingale property. Specifically, we con- 
struct, analyse and classify three new families of nonlinear diffusion models with affine 
drift that arise from the squared Bessel process (the Bessel family), the C1R process 
(the confluent hypergeometric family), and the Ornstein-Uhlenbeck diffusion (the OU 
family) . 

Introduction 

A solvable continuous-time stochastic process can be denned as a process for which the 
transition probability density function is obtainable in analytically closed-form. This kind 
of solvability permits us to precisely simulate paths of the process from its exact sample 
distribution. Another important benefit is that some standard methods such as spectral 
expansions can be applied to obtain closed-form analytical formulas for the transition den- 
sities of regular as well as killed solvable diffusions. As a result of that, one can derive other 
fundamental quantities such as first-hitting time distributions as is done in |13j and [5]. 

The set of diffusion processes that are, on the one hand, tractable and applicable for 
mathematical modelling and, on the other hand, exactly solvable is not so vast. This set 
includes mostly linear diffusion processes or diffusions with power nonlinearity (see [4] and 
P2] for a comprehensive review of such diffusions). There are two main tools that allow 
us to construct new solvable diffusion processes. The first is related to a measure change 
on a chosen underlying diffusion and the second involves a change of variable or smooth 
monotonic mapping (the Ito formula). In recent years, a new approach that combines special 
measure changes with nonlinear variable transformations was introduced for uncovering 
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various new families of exactly solvable driftless diffusion models [2J |3J [BJ [TT] . The method 
has been coined as "diffusion canonical transformation", partly in light of the fact that a 
diffusion process is reduced, or mapped, onto another underlying (simpler) diffusion process. 
In Subsection |1.1| we state the standard Green's function methodology and provide the 
basic background formalism and general assumptions that define the underlying diffusions. 
Subsection |I . 2| presents the diffusion canonical transformation method. Detailed discussions 
of other canonical aspects of the transformation methodology can be found in [5]. 

In this paper we discover that a relatively simple change in the nonlinear variable trans- 
formation allows us to construct a whole new class of diffusions (F t )t>o with linear drift and 
nonlinear diffusion coefficient functions defined by the time-homogeneous stochastic differ- 
ential equation (SDE) of the form dF t = (a + bF t )dt + <r{F t )dWt- As a result of our study, 
we upgrade three main families of time-homogeneous diffusions named as Bessel, conflu- 
ent hypergeometric, and Ornstein-Uhlenbeck families. These processes, which we also refer 
to as F-diffusions, arise via the diffusion canonical transformation method by respectively 
choosing the standard squared Bessel process, CIR process, and Ornstein-Uhlenbeck process 
as underlying diffusions (see Section [3]) . It is important that these three families include 
and recover all the corresponding families of driftless diffusions obtained previously in [BJ as 
special subfamilies. Moreover, the new affine models inherit some of the important salient 
properties of their driftless counterparts. 

One immediate application of such diffusions is asset pricing in finance (when a = 
and b is a constant interest rate). In [BJ, we show that these three families generate local 
volatility profiles with varied pronounced smiles and skews. Moreover, we prove that there 
exists a risk-neutral measure so that the driftless asset price process is a martingale. Option 
pricing under these models is straightforward and is well studied in J5j0[8]. Clearly, pricing 
of European options is reduced to the evaluation of a definite integral. Since the first- 
hitting time distributions for killed solvable diffusions is available in the form of a spectral 
expansion, no-arbitrage prices of barrier and lookback options are available in closed form 
as well. As well, discretely-monitored path-dependent options can be evaluated by using a 
path integral approach. In doing, so we use one important observation that the distribution 
of a bridge asset price process is reduced to the distribution of the corresponding underlying 
bridge diffusion (e.g. the Bessel, CIR, or Ornstein-Uhlenbeck bridge processes). Moreover, 
subfamilies of diffusions belonging to the Bessel and confluent hypergeometric families admit 
absorption at a boundary point, so they can naturally be used in derivatives pricing under 
credit risk. Last, but not least, the newly developed F-diffusions are also applicable to 
interest rate modelling. Although the diffusions with affine drift a + bF considered here are 
not mean-reverting, they still can be useful as short-term interest rate models thanks to the 
nonlinear local volatility function. Alternatively, by using a suitable change of variable, we 
are able to obtain models with nonlinear mean-reversion that permit closed-form solutions 
for valuing zero-coupon bonds. 

One of the main contributions of this paper (in comparison with [6j and references 
therein) is the analysis and classification of allowable mapping functions that lead to new 
diffusions with affine drift. Subsection |1.3| presents some results for proving the smooth 
monotonicity property of mappings that lead to F-diffusions with affine drift where 6^0. 
Another main result, presented in Subsection |1.4| is the complete boundary classification of 
such newly constructed regular diffusions. Finally, we analyse whether or not such diffusions 
{F t )t>o "preserve" the drift rate, i.e. whether ^E[F T | F t ] = E[a + bF T \F t ],t<T, holds. 
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This property can be viewed as a generalization of the martingale property for driftless 
processes, where E[F T | F t ] = F u < t < T, and, therefore, ^E[F T | F t ] = holds. 
Therefore, for the special case with a — we prove that the discounted asset price process 
(e~ bt F t is a martingale. With the aim to provide a rigorous framework for the above 
characterizations of the families of F-diffusions, in Section [2] we present a theorem that 
state easy-to-implement limit (asymptotic) conditions for verifying the conservation of the 
expectation rate of the ^-diffusions. 

In Section [4] we apply the theory of Section [l] and Section [2] to the above-named three 
chosen underlying solvable diffusions defined in Section|3] In doing so, we give a construction, 
probabilistic classification and analysis of the resulting new three families of ^-diffusion 
models. In Section [472] we discuss all the possible monotonic maps that lead to nonlinear 
^-diffusions with affinc drift a + bF where b ^ 0. Some conclusions are drawn in Section [5] 
The Appendix contains useful asymptotic properties of certain fundamental solutions that 
characterise the three known underlying diffusions. 

1 General Construction of Nonlinear Solvable Diffu- 
sions with Affine Drift 

1.1 Underlying diffusion 

Let (X t ) t >Q be a one-dimensional time-homogeneous regular diffusion on X = (I, r), -co < 
I < r < oo, defined by its infinitesimal generator: 

g x f = (gf)( x )^^ 2 (x)f"(x) + \(x)f(x), xei. (i.i) 

The functions X(x) and v(x) denote, respectively, the (infinitesimal) drift and diffusion 
coefficients of the process. Throughout we assume that the functions X(x), X'(x), v{x) and 
v"{x) are continuous on the open interval I and that v(x) is strictly positive on I. 

Let px(t;xo,x) denote a transition probability density function (PDF) for the diffusion 
specified by the conditional probability that X t is contained in an interval DCI: ¥{X s+t € 
D\X S = .To} = J D Px{t; xq, x) dx, \/xq £ X, t, s > 0. Since the drift and diffusion coefficients 
are assumed to have no explicit time dependence, then p = px(t; xo,x) satisfies the time- 
homogeneous forward Kolmogorov equation ^ = C x p and backward Kolmogorov equation 
§f — Qx P subject to the Dirac delta function initial condition p(0+; £o, £) — S(x — xq). 
Here C x denotes the Fokker-Planck differential operator (see, e.g., [10]). 

The diffusion (X t )t>o can also be defined by its speed measure M(dx) and scale func- 
tion S(x) (see, e.g., [1]). For diffusions considered here, these characteristics are absolutely 
continuous with respect to the Lebesgue measure and have smooth derivatives. The speed 
and scale density functions m(a;) = M ^ x x ^ and s{x) — dS } x ^ are as follows: 

/ f x 2\(z) \ , . 2 . . 

six) = exp - / dz and m(a;) = „.,., . (1.2) 

V J v \ z ) J v 2 (x)s(x) 

The functions s(x) and m(x) are hence continuous and strictly positive on X. The infinites- 
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imal generator and Fokker-Planck differential operator can be re-written in compact form: 



m(x) \ s(x) J \ s ( x ) \ m (' T ) 

As is well known, there are two linearly independent fundamental solutions (pf and <pj 
of the ordinary differential equation (ODE) 

(gcp)(x) = sp(x),seC,xel, (1.3) 

such that for positive real values s = p > the solutions ipp~(x) and <fp(x) are correspond- 
ingly increasing and decreasing functions of x (see, e.g., @]). The functions ip^(x), p > 0, 
are convex and finite in the interior of the domain T. 

In this paper we consider underlying diffusions (X t ) t > for which the endpoints I and r 
are either unattainable or regular instantaneously reflecting boundaries (see Feller's classi- 
fication in nun]). 

Let Sx,S2 € C with positive real parts, then the functions ip^ satisfy the asymptotic 
relations 

ft, ( x ) ft, ( x ) 
lim — ^ = and lim — ^ = oo. (1-4) 

x^l+ Lp S2 (x) x->r- (p a3 (x) 

Moreover, the square integrability conditions 

(<pt>Vp)Vrt < 00 and (Vp'^p W) < oo (1.5) 

hold for p > and i£l. Throughout this paper we conveniently define the inner product 
of two functions /, g with m on a closed interval [a, 6] as 

{f,9)[a,b]= f(x)g(x)m(x)dx (1.6) 

J a 

and (f,g)(a,b] = lim (f, g) [e . b] , (f,g) [a , b ) = lim {f,g)[ a ,e\, WfWLb) ~ (/>/)(«,&)• 
The Wronskian of the fundamental solutions can be computed as follows: 

W[<p-,<p+]{x) 4 ^(z)^j^ - = ".«(*) , (1.7) 

where id s is a constant w.r.t. a; and u> p > for real p > 0. 

The time-independent Green's function G(x,Xo,s) is defined by the Laplace transform 
of px(t; xq, x) with respect to time t. Consequently, a solution for the PDF px(t; xq, x) 
follows by Laplace inversion, as given by the Bromwich integral. For an in-depth discussion 
of this inversion and how it leads to various spectral expansions for the transition PDF, see 
[H] and [5]. 

The Green's function G, for x,x$ £ X, is written in terms of the functions (pf, <pj, and 
m(x) in the standard form jl]: 

G{x, x , s) = w~ 1 m(x)ip+ (x^ip- (&>), (1.8) 

where x < = min{a;, xo} and x> = max{a;, io}- The boundary conditions for G follow from 
those satisfied by the appropriate solutions ipj and tpj. 
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1.2 Diffusion Canonical Transformation 

Consider a class of one-dimensional time-homogeneous regular diffusions (X^ [ )t>o £ 1 with 
infinitesimal generator 

1 / v! (x) \ 

{Q (P) f){x) ± ^ 2 {x)f\x) + X(x) + v 2 {x)^\ f(x) , (f .9) 

2 V u p (x)J 



where X(x) and v(x) are the infinitesimal drift and diffusion coefficients of the underlying 
process A t defined by (jXTTJ) . A strictly positive generating function u = u p {x) is a linear 



combination of the fundamental solutions ip^ 



u p (x) = qxipp(x) +q 2 ip p (x), (f.fO) 

with parameters qi,q2>0 and at least one of them being strictly positive. Throughout, p is 
a real positive parameter. The speed and scale densities for an Jf ^-diffusion are given in 
terms of those for the underlying A-diffusion: 

s( x ] 

m p (x) = u 2 (x)m(x) and s p (x) = „„' , . (1-H) 

v£{x) 



By comparing the generators (1.1) and (f.9), observe that A( p )-diffusions can also be 
viewed as arising from the underlying A-diffusion by the application of a measure change. 
In fact, the -diffusion can be realized from the X-diffusion upon employing a so-called 
Doob-ft, transform, where h = u p , or p-excessive transform. For more details on this and 
its connection to Girsanov's Theorem, see [5]. Both processes are regular on the same state 
space X = (l,r). A transition density Px\t; xq, x) for an A^-diffusion is then related to 
that for the A-diffusion (in the same original P-measure) by 

p { p(t;x ,x) = e- pt ^\px(t;x ,x), x,x o el,t>0. (1.12) 

Up{X ) 



Its important to note that this relation (as well as equation ( 1.15 1 below) also holds for killed 
diffusions. In such cases, regular killing boundary conditions are (additionally) imposed at 
one or two interior points and the transition PDFs for the A and A^ killed diffusions 
are defined on the same subspace of T (see [5] for theory and applications of various killed 
diffusions) . 

We now consider ^-diffusions {Ft = F(X^),t > 0} defined by strictly monotonic real- 
valued mapping F — F(x) with F', F" continuous on X and having infinitesimal generator 

{Qvh)(F) = i a 2 (F)h" (F) + (a + bF)h'(F), F e If = (F« , fW) . (1.13) 

(Ft)t>o is a regular diffusion on If with endpoints F^ — min{F(i+), F(r— )} and F^ = 
max{F(Z+), F(r— )}. The diffusion coefficient function <r(F), as given by (1.23) below, is 
generally nonlinear and where a and b are arbitrary real constants such that 6 = implies 
a = 0. The corresponding scale and speed densities are 
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Several families of analytically solvable state dependent volatility models can be derived 
from known underlying diffusion processes defined by ( 1.1 1. We refer to this construction as 
the "diffusion canonical transformation" methodology. For driftless diffusions, see [21 El E]- 



Extensions and applications to new families with linear drift (a = 0. b 0) are discussed in 

The transition PDF pp for an _F-diffusion (F t ) t >o is related to the transition PDF for 
the underlying X (or X^) diffusion as follows: 



PF(t;F ,F) 



^^;X(F ),X(F)) 

KX(F)) u p (X(F)) _ pt 
a(F) u p (X(F )) 



(1.15) 



Px (t;X(F ),X(F)) . 



where F, F £ lp, t > 0. Here X = F 1 is the inverse map so that 

u(X(F)) 



\X'(F)\ = 



a(F) 



This methodology was originally developed for driftless F-diffusions where a = b = 0. 
For such cases, the volatility function has the form 



= roKsMs) x 
u 2 p (x) 



X(F) , ao > . 



(1.16) 



The map F(x) (that solves equation (1.181 below for the special case a = b = 0) admits the 
general quotient form: 

c 1 (p+(x)+c 2 (p~(x) 



F(x) 



(1.17) 



Qlft(x) +q2<p p (x) ' 

where Ci, C2,<fi, <fe € E. are parameters such that q\Ci — q2^\ ^ 0. Several families of F- 
diffusions arising from various choices of underlying diffusions (such as the squared Bessel, 
Ornstein-Uhlenbeck, CIR and Jacobi processes) are studied in [21 [6l [5] . In this paper, we 
expand the families constructed by finding new mapping functions F so that the newly 
obtained diffusions have generator in (1.131, i.e. they satisfy an SDE with affine drift: 

dF t = (a + bF t )dt + a{F t )dW t ■ 

Let F be a strictly monotonic twice continuously differentiable function. Applying the 
Ito formula, the process F t = F(X { t p) ) satisfies the SDE 

dF t - {G (p) F){X{Ft))dt+v{X{FtW{X{F t ))\dWt . 

Clearly, by equating these two SDEs, the mapping function F(a;) leads to a process with 
affinc (linear) drift defined by the generator in (1.13) and solves the following 2nd order 
linear nonhomogeneous ODE: 

(G ip) f)(x) = a + bF(x). (1.18) 
The strictly positive diffusion coefficient is then o{F) — v(X(F))\F'(X(F))\, or 

o-(F(x)) = v(x)\F'(x)\. (1.19) 
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Lemma 1. Let 6^0 and p, p + b > hold. Then the solution to equation (1.18) takes the 
general form 

a c lV >+ +b (x) + c 2 ip- +b (x) _ a v p+b {x) 

b u p (x) 



F(x) = - 



(1.20) 



b ' qitpp~(x) +q 2 <p P (x) 
where C\ and C2 are arbitrary real constants. 

Proof. The numerator function v = v p+b {x), denned in (1.201, is a linear combination of 
and hence solves Qv = (p + b)ii. The positive function u — u p {x) in the denom- 
inator solves Qu — pit. Now differentiating and using the identity G (p) f(x) = Qf{x) + 
v 2 (x)(u' (x)/u p (x))f'(x), we readily have 



?0) 



V 


1 












Gv- 






u 


u 




u 1 





- [(p + b)v- ~{pu) 
u 



bl 
u 



Hence, 5 is a general solution to the corresponding homogeneous ODE (eq. ( 1.18 1 for a = 0) 

W[<pi +b ,<p:+ b K*) 



since the Wronskian W 



(x) = 



"p+b^ P +b 

u 2 (x) 



function F p (x) = —a/b is a particular solution of (1.18). 



7^ from \1.7\. The constant 

□ 



The derivative of the mapping in (|1.20|) is simply 
F'(x) = 



v' p+h {x)u p {x) - u' p {x)v p+b {x) 



W{x) 



(1.21) 



u 2 p (x) u 2 p {x) 

where we define the Wronskian 

W(x) = W(x;p,p + b) 4 W[u p ,v p+b ](x) . (1.22) 

Assuming F is strictly monotonic, and given an X-diffusion, the F-diffusion coefficient func- 
tion a(F) is then given by substituting (1.211 into (1.191: 

v(x)\W(x)\ 



a(F) 



u 2 p (x) 



X(F) , Fel F . 



(1.23) 



We note that this expression holds for all parameter choices a, b except when a ^ and 
6 = 0. For the latter special case (i.e. with constant nonzero drift function) equation (1.18) 
reads G^ F(x) = a and hence simply reduces to a linear first order ODE in F'. Solving leads 
to various monotonic maps which in turn give rise to nonzero constant drift F-diffusions 
with various nonlinear specifications for the diffusion coefficients. In this paper, we shall not 
discuss the details of such special families as we focus on linear drift functions with 6^0. 



1.3 Monotonic Maps 

The map F : X — * Xp in ( |1.20[ ) does not generally satisfy F'(x) 0. To guarantee that 
(Ft)t>o is a regular diffusion on If = (F^, f! r !) the map F has to be strictly monotonic. 
Then F'(x) 7^ and hence the diffusion coefficient function a(F) is strictly positive on Xp. 
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From ( 1.21 ) we observe that the sign of F' equals the sign of W(x). Using the represen- 



tations of Up and v p +b in terms of ipp , f p+b gives 

W(x) = q 1 c 1 W[^+^ + p+h ]{x)+q 1 c 2 W[^+, V - p+h ]{x) 

+ q 2 c Y W[ip~, <pp~ +b ](x) + q 2 c 2 W[ip~ , <p~ +b ]{x) 



(1.24) 



There are two important cases where F is strictly monotonic. Recall that the fundamental 
solutions tp~^ (x) and ip~{x) are correspondingly strictly increasing and decreasing functions 



of x. Therefore, the ratios 



and fI±M 



are strictly increasing and decreasing func- 



tions, respectively. In particular, we have the strict inequalities W[ipp~ ,ip p+b \{x) < and 
W[ip~ , Vp+bK^) > 0. Thus, the two choices of parameters c 2 = q\ = 0, c\jq 2 = ±c or 
c\ = q 2 = 0, c 2 /qi = ±c in equation (1.20) lead to dual subfamilies of strictly monotonic 
maps defined by F = F^: 

+ 



fW(x) 



7+fC 



<Pp+b( x ) 



(1.25) 



where e = ±1 and c > is constant. Other parameter choices that lead to other families 
of monotonic maps are discussed in Section |4.2| The following propositions are useful in 
verifying whether or not W(x), and hence f'(x), changes sign in X. 



Proposition 1. The Wronskian W(x) in 1 1.22) satisfies the first order linear ODE l>v 2 (x)W' (x)- 
\{x)W(x) — bupVp+b! and for any x, Xq G T the solution admits the following representation: 



W{x) _ W(x ) 

S(x) S(x ) 



b / m(y)up(y)v p+b (y)dy . 



(1.26) 



Proof. The proof follows by direct verification upon using Qu p = pu p and Qv p+ b — {p + 
b)v p+b . □ 



(1.71 and (1.241 



b = 0, the function V ^}- x ) is constant and from 



Note that for the driftless case, with a 

W(x) = (ciq 2 — c 2 qi)w p s(x). Therefore, equation (j 1 . 1 6|) defining a(F) 



for all families of driftless _F-diffusions is recovered as a particular case of the more general 



specification given by equation (1.23). Moreover, the specification in (1.23) gives rise to a 



class of state dependent volatility models that generally can also have a dependence on the 
drift parameters a and b. 



Proposition 2. Assume that c\ and c 2 in $1.20 ) are both nonnegative or nonpositive (with 
at least one of them being nonzero)and that W(x) in ( 1 . 22\ i preserves its sign as x approaches 
either endpoint I or r; that is, sign(M / (Z+)) = sign (W(r— )) . Then Vi e I W(x) ^ 0, i.e. 
W(x) is either strictly positive or negative on T. 

Proof. Under the assumed conditions on c\ and c 2 , the function Wp+f, is either strictly positive 
or strictly negative on X. Recall that u p is strictly positive. Hence, the function W (x) / s(x) 
given by ( |1 .26 1 is monotonic in x and, since s(x) > 0, W(x) has at most one zero in T. Then 
sign(W '(l+)Y~= sign(l4 7 (r-)) implies W(x) ^ 0, i.e. either W(x) > or W(x) < for all 

x el. □ 
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Proposition 3. Assume that c\ and c 2 for the map F in (1.20) are both nonzero and have 
opposite signs. 

(i) lfb>0, then W(x) is either strictly positive or negative on X. 

(ii) Letb<0 and xeX. If sign(W (1+)) = sign(W(r-)) = +1 and c x > 0, then W(x) > 0. 
Ifsign(W(l+)) = siga(W(r-)) = -1 and c x < 0, then W(x) < 0. 

Proof. By definition of v p+ b, and the fact that f^ +b and <p p+b are respectively increasing and 
decreasing functions, we have v' p+b (x) > (or i' p+b (x) < 0) when c\ > 0, c 2 < (or cj < 0, 
c 2 > 0), i.e. Vp+b{x) is either a strictly increasing (or decreasing) function. Moreover, if 



follows from such monotonicity and the properties in (1.4 1 that v p+ b(x) has exactly one 
zero, at x = x € X where ip+ +b (xo) / <p~ +b (x ) = \c 2 /ci\. Then W( x ) = u p {x )v' p+b {x Q ) is 



accordingly strictly positive (or negative). Setting xq = xq in (1.261 now gives 

- ' 1 m(y)u p (y)\v p+b (y)\dy 



S(a;) s{xo) Jmin(x ,x 

where e = +1(-1) if c x > (c x < 0). Hence, if b > then either W(x) > or W(x) < 0, 
for all x e X, in the respective cases. If b < 0, then W(x), as given in the last expression, 
can either have no zeros or at most two zeroes in X. W(x) has no zeros if and only if 
sign(T / K(io)) = sign(W / (Z+)) = sign(H / (r— )). This hence proves part (ii) for the respective 
cases. □ 

As follows from the above propositions, to determine the monotonicity of a map F, of 
interest is the asymptotic behaviour of W{x), as x approaches endpoint I or r. Below 
we consider three families of F-diffusions arising from three separate underlying diffusions; 
namely, the squared Bessel (SQB) process, the Cox-Ingersoll-Ross (CIR) model, and the 
Ornstein-Uhlenbeck (OU) process. For all these families, the asymptotic properties of the 
fundamental solutions ipf and of the corresponding Wronskian functions are presented in 
Appendix [X] 

1.4 Boundary Classification for F-Diffusions 

Given an underly ing A-diffusion and p > 0, any regular diffusion (X t )t>o € (I, r) with 



generator g {p) in KM falls into one of three general families: 

(i) {qi = 0,q 2 > 0} where u p (x) = q 2 <p p (x), 

(ii) {qi > 0,q 2 = 0} where u p (x) = qi(p+(x), 

(iii) {qi,q 2 > 0} where u p (x) = qi^(x) + q 2 ip p (x). 

Lemma 2. The above three fam ilies (i)-(i ii) of regular diffusions (x) p ') t >o on {l,r) with 



generator Q( p \ p > 0, defined by (1.9) and (1.10) have the following boundary classification: 

(i) q\ = 0,q 2 > 0: The boundary I is attracting natural if (<pp,tPp)(i tX ] — °°; e2: ^ */ 
(<P P ~,<Pp)(l,x] < oo and(<p p ,<p p )(i, x ] = oo, and is otherwise regular when (ip p , V~)(i, x ] < 
oo . The boundary r is non-attracting natural if (tpj , ip~ )[ x , r ) = 00 an d * s otherwise 
entrance. 
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(ii) q\ > 0,(72 = 0: The boundary r is attracting natural if (ip~£ , tp~ )\ x , r ) — °°> ex it if 
(tpt, ip~p~)[x,r) < 00 an d { l Pp> ft )[x,r) — °°; an d * s otherwise regular when ((ft, <pp~ )[ x ,r) < 
00. The boundary I is non- attracting natural if (ipti <Pp~)(l,x] = 00 ar id is otherwise en- 
trance. 

(Hi) q\ > 0, qi > 0: The boundary I has the same classification as in (i) and r has the same 
classification as in (ii). 

Proof. See 0. □ 

According to this Lemma, the left and right boundary classification follows simply from 
the asymptotics of m(x)<Pp(x)ip~(x) as x — ► 1+ and x — > 1 — , m(x)[ip+(x)} 2 as x — > r— and 
m(x)[ip~ (x)} 2 as x — > i+. 

Theorem 1. XTie boundary classification for an F -diffusion is equivalent to the correspond- 
ing jW -diffusion. 



Proof. This follows by the diffeomorphism x[ p ^ — > F f = F(X( P- '). See [§] for details 



□ 



2 Drift Rate Conservation and Martingale Property for 
F-Diffusions 



For any time-homogeneous _F-diffusion defined by ( 1.13 1, we introduce the rate of change of 
the conditional expectation: 



d_ 

at 



E[F T+t I F T = Y] = [ F^-p F (t;Y,F)dF, Y e If, * > 0, r > . 



The PDF pi? given by (1.15) satisfies the forward Kolmogorov equation 

dpF d ( 1 d ( pf 
~cW = dF [s^dF W 



(2.1) 



(2.2) 



with densities in (1.14) now defined in terms of those for the -diffusion: 



m F (F) = m p (X(F)) \X'(F)\ and s F (F) = s p (X(F)) \X'(F)\ . (2.3) 
Using (2.2 1 within (2.1|, and applying integration by parts together with the derivative 



d 

(IF 



(s^Fyj = ( a + bF)mF(F), gives the rate in (2.1l expressed as a sum: 




r+t 



F T = Y}= / (a + bF) PF (t; Y, F) dF + £(Y,t) . 
Jfw 



dt Jf i, 
The "bias" term £(Y,t) is given by the difference of two limits: 



(2.4) 



£(Y,t) 



F f PF (t;Y,F) 



1 p F (t;Y,F) 



s F (F) dF V m F (F) J s F (F) m F (F) 
F(x) d (p ( Z\t;y,x)\ F'(x) p ( ?\t; y, x) 



n F=F^ 



F=F(0 



s p (x) dx \ m„(a;) 



s P {x) m p (x) 



(2.5) 



x=l+ 
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The last expression follows by changing variables x — X(F),y = X(Y") and by combining 
( [23| and ( [Lis) . 

Consider the case where £(Y,t) = 0, for any F T = Y E Zp,£ > 0. Then, for such 



diffusions one obtains a simple representation for the rate of change in (2.4 1 as 



0_ 

Ft 



E[F T+t | F T ] = aE[t FT+teIf | F T ] + bD[F T+t \ F T ] . 



(2.6) 



i.e. the rate of change of the conditional expectation equals the conditional expectation of 
the drift function for the process. If (F t )t>o conserves probability i.e. if ~E[1p t£ x F I F r ] = 
F{F T+t € If\F T eI F } = l for all r > 0, t > 0, or if a = holds, then E(t) = E[F r+t | F T = 
Y] satisfies the trivial linear ODE E'(t) = a + bE(t), subject to E(0) = Y. Let a = 0, then 
there exists a well-known geometric drift solution E(t) — Ye bt , i.e. 'E\F T+t \ F T = Y] = 
e ht Y for all Y G 2f,t > 0,t > 0. In other words, the "discounted process" (e~ bt F t )t>o is 
a martingale in case a = and £ = 0. Note also that the discounted process is a strict 
supermartingale (submartingale) when a — and £ < (£ > 0) for all £ > 0. Setting r = 
recovers the unconditional expectation E[F t ] = i 7 be fct . 



By Laplace transforming (2.5 1 with respect to time £, changing the order of integration 



and differentiation, while employing the assumptions formulated in Subsection |1.1[ and 
using (1.8 1, (1.121 and (1.211, the following theorem is readily proven. A similar proof for 



driftlcss F-diffusions is given in [B]. 



Theorem 2. Let F be a monotonic map given by (1.20) for 6^0 an d by (1.17) for 

a = b = 0. The diffusion F t — F(X^),t > 0, with (X^) t >Q defined by (1.9), conserves 
the expectation rate, 
conditions hold: 



lim 



lim 



i.e. the relation (2.6) is true if and only if the following boundary 
W[up,tp+ +S ](x) W[u p ,v p+b ](x) <fp +s (x)~ 



F(x) 



F(x) 



s(x) 



s(x) 



u p {x) 



W[up,tp p+S ](x) W[u p ,i p+b ](x) <p p+s (x) 



= 0, 



= 0, 



s{x) s{x) u p {x) 

for all complex-valued s such that Re s > c, where c is some real constant. 



(2.7) 



(2.8) 



3 Three Choices of Underlying Solvable Diffusions 

3.1 The Squared Bessel Process 

We specifically consider an underlying Ao-dimensional squared Bessel (SQB) process obeying 
the SDE dX t = X n dt + v \fX~ t dW t with constants v > and A > 0. The scale and speed 
densities for the diffusion are six) — x^^ 1 and m(x) — \x^ , where u = ^ — 1. Here we 

take Ao > hence [i > 0, where the SQB process is defined on X — (0,oo). The point 

x = is an entrance and x — oo is an attracting natural boundary. For the regular diffusion 
on X the transition density function is 
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The generating function u p {x) is a linear combination of 

ip+(x)=x-^ 2 I^2y/2px/u ) and <pj(x) = x~^ 2 K^2^2px/v ) , (3.2) 
where Ift(z) and K^(z) are the modified Bessel functions (of order jj) of the first and second 



kind, respectively (for definitions and properties see [T]). The pair (pf(x) satisfies (1.7) 
where w s — 1/2. For Si,s 2 £ C, all Wronskians W[iff i: (pf 2 ](x) and W[ip~ i ,ip~^ 2 ]{x) are 
readily obtained using differential recurrences zl'^z) = ^^(z) + zl fl+ i(z) and zK'^z) — 
nK,t(z) - zKn + i(z). 

3.2 The CIR Process 

Consider the Cox-Ingerssol-Ross (CIR) process [9] (X t ) t >o £ X = (0, oo) solving the SDE 
dX t — (Ao — \\X t )dt + v^\j 'X t dWt, where A > 0, Ai > 0, and v > 0. Here we take 
= 2^ — 1 > 0. The endpoint x = oo is non-attracting natural and the origin is an 

entrance boundary. The respective scale and speed densities are s{x) — x~^ e KX and 
mix) = -^x^e^^, where k = > 0. For the regular diffusion on X the transition PDF is 



(3.3) 



p x (t;x ,x)=c t e^ t ( ^— J e -<=*(-« Xlt +*o) ^ ( 2ci ^) , 

where ct = n/(e Xlt — 1). 

The generating function u p (x) is a linear combination of the functions 

<Pp(x) = M.(v,/j,+ 1,kx) and tp~ (x) = U (v, p + 1, kx) (3-4) 

where v = p/\\ > 0. The functions M(a, 6, z) and U(a,b,z) are confluent hypergeometric 



functions, i.e. the standard Kummer functions (see pQ). The pair ff(x) satisfies (1.7) 
where w s = For si,s 2 € C, W[^ 1 ,^](a;) and W[<pJ ± , <pf 2 ](x) are obtained 

using differential recurrences 4-M.(a,b,z) — (a/b)M(a + 1,6 + l,z) and j-U(a,b,z) = 
-aU{a + 1,6+1, z). 

3.3 The Ornstein-Uhlenbeck Process 

The Ornstein-Uhlenbeck (OU) process solves the SDE dX t = (Ao — XiX t )dt + i/odWt, where 
Ao,Ai > 0. Both boundaries, I — — oo and r = oo, of the state space X = (—00,00) 
are non-attracting natural for all choices of parameters. Without loss in generality we set 
Ao = 0. Otherwise we can consider the shifted process Y t = X t ~ ^ and the formulas 
follow by simply shifting x — > x — y 2 -, xq xq — ^+ The speed and scale densities are 



s(x) — e KX I 2 and m{x) — (2/v 2 )e KX ~/ 2 . As in Subsection 3.2 we define the positive 



= 2 ^r- , v = #- . The regular OU diffusion on X has the transition PDF 

Kg ' Al & 



= \/ 2tt(1 -l-^i') ° XP (" "2(1 -l-Kit) ' ' 
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The generating function u p {x) is a linear combination of 

<p~(x) =e Kx2 / i D_ v (^Kx) and <p+{x) = <p~(-x) (3.6) 

where D_ v (x) is Whittaker's parabolic cylin der f unction (see £Q for definitions and prop- 
erties). The Wronskian constant in equation ( jl ,7| > is given by w s — r ffi^ ) ■ For si,S2 G C, 

W[y>^ , <pf 2 ] (x) and W[y>7 , <pf 2 ]{x) are obtained using differential recurrences 4^D- V {z) = 
-{z/2)D- v (z) -vD-v-x \z). 



4 The Bessel, Confluent Hypergeometric, and Ornstein- 
Uhlenbeck Families of Affine Diffusions 

4.1 Three Main Families of Affine Drift Diffusions: Classification 
and Properties 

Using the general construction presented in Section [T] and the three underlying X-diffusions 
from Secti on [3] we can construct three new families of F-diffusions with affine drift as de- 
fined by (1.13) and (1.231. For all such new families, the transition PDF pp is given in 
analytically closed form for regular as well as killed _F-diffusions. This partly explains why 
these processes are also termed "solvable" diffusions. Generally, a transition PDF pp is 



given by (1.15), in terms of the corresponding transition PDF px of the underlying X- 



diffusion, where the monotonic map F(x) (with its inverse map X(F)) is given in (1.201 and 
the generating function u p is given in (1.10 1. By choosing the SQB, CIR, or OU diffusion 
as underlying X-diffusion, we respectively obtain the Bessel, confluent hypergeometric, or 
Ornstein-Uhlenbeck families of F -diffusions with affine drift. Some properties and classifi- 
cation of the corresponding driftless diffusions are discussed in [BJ. For regular diffusions 
on the respective state spaces 1 = (0, oo) or X = (— oo, oo), the transition PDF px is given 
by equation (3.1|, (3.3), or (3.51, respectively. For corresponding diffusions with imposed 



killing at one or two interior points of the state space, the transition densities px are given 
in terms of closed- form spectral expansions and these, in turn, lead to closed- form spectral 
expansions for transition PDFs pp, as well as other fundamental quantities such as first hit- 
ting (passage) time densities for i^-diffusions (see jS]). The fundamental solutions ^ used 



in the map and generating function are given by either (3.2), (3.4 1, or (3.6 1, accordingly. 



Each i^-diffusion is described by the set of parameters that are inherited from the chosen 
underlying diffusion. In addition to that set, the nonnegative parameters q±, and p > 
are added thanks to the measure change X — ► X^ p \ Two parameters a and b describe the 
affine drift coefficient in (1.131. Finally, up to two other parameters c\ and C2 are used 
in the map function. As observed from the diffusion coefficient function o~(F) in (1.231, 
the combination of all such parameters make the new diffusions quite flexible for modelling 



various stochastic processes. Different choices of monotonic maps F = F± , i = 1, . . . , 5, lead 
to different _F-diffusions. In any case, the diffusion function is specified by equation (1.231. 

The choice 
finance 



± leads to _F-diffusions that have applications in asset pricing in 
In particular, the dual maps defined by (1.251 with e = +1 give rise to sets of dual 
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subfamilies (i) and (ii) of affine drift .F-diffusions with respective volatility specification: 



a(F) = cv{x) 



Wp{x)Y 

w[y>; +b ,y>+Kx) 



(ii) 



= cv(x) 



Wp(xW 

Vp+bi^^ft^ 



<Pp~(x) 
<Pp( x ) 



(i) 
(ii) 



(4.1) 



where x — X(F) = F is the unique inverse map for the respective subfamilies (i) F(x) 

. c - : 'iii/i riii i— r ' ' i ... / 



Vp ( x ) 



and (ii) F(:r) 



a , fp+b\ 



Both subfamilies have regular state space 



F € (— f,co). Computing the Wronskians in equation (4.1l gives the diffusion coefficient 
function for three dual subfamilies as follows: 



a{F) = cV2 



K,(^V2px-) 

'±y/2( P + b)x) 



-bKu 



(1) 

(ii) 



(4.2) 



for the Besscl family, 



cr(F) = ckvo\/x 



v M[y+ ^ ,h+1,kx)U(v+1,h+2,kx) (v+ M (v+ ^-+1,^+2, «i) 



hi 2 (w,/i+l,Ka:) 



U(v,p+1,KX) 



vu(v+-^,p,+l,Kx^M(l+v,fi+2,Kx) (v+-^)u(v+^-+1,ij,+2,kx^ 



O+l) M?(v,h+1,kx) 



A4(v,p+l,Kx) 



(i) 
(ii) 



for the confluent hypergeometric family, and 

( (p+b) D_ (v + j^ ) + 1 {-^kx) pD_ (v+ b ) (-\/Kx)D_( ll+1) ( s /Kx) 

~2~ 



a(F) = c 



D- v (^kx) ~ D 2 _ v {^/kx) 

(p+b) D_ (tj+jL+1) ( v ^a;) p D_ (t)+ 6 } {s/k x)D_( v+1 ) (-y/Ux) 



D_„(-V^a:) 



(i) 
(ii) 



(4.3) 



(4.4) 



for the OU family. In all cases, x — X(F) = F~ 1 (F) using the respective map for (i) or 
(ii). Subsets of these families with a = 0, b ^ have regular state space F € (0, oo) and are 
useful for modelling asset prices in finance. Figures [l] and [2] display some computed curves 
of the local volatility function ai oc (F) :— a(F)/F for the subfamilies in equations (4.2 i), 
(4.3 i) and (4.4) when a — 0,6 ^ 0. We note that the dual OU subfamilies (i) and (ii) 
coalesce when a = 0. As seen in Figures [T] and [2j by adjusting parameters, the models are 
all readily calibrated to attain a prescribed level of local volatility for a given value of F. 
The respective sets of freely adjustable positive parameters for the dual Bessel, Confluent 
and OU subfamilies are: (p, 6, p,, c, vq), (p, 6, k, /i, c, vq) and (p, 6, k, c, vq). 

The following two theorems describe the properties of the Bessel, confluent hypergeo- 
metric, and Ornstein-Uhlcnbeck families of F-diffusions with affine drift. These theorems 
generalize the results obtained previously for the driftless case (see [B]). 

Theorem 3. The Bessel, Confluent and OU families of regular diffusions F t — F(X^), 
with affine drift and monotonic map defined as in Theorem^ have the following boundary 
classification. 
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Bessel and Confluent families: The endpoint F = F(0+) is entrance if q 2 — 0; is 
exit if q 2 > and \x > f ; is a regular killing boundary if q 2 > and < fi < 1. The 
endpoint F = F(oo) is non- attracting natural if qi = 0; is attracting natural if qi > 0. 

OU family: The endpoint F — F(— oo) is non-attracting natural if q2 — and is 
attracting natural if q 2 > 0. The endpoint F = F(oo) is non-attracting natural if 
qi = and is attracting natural if qi > 0. 

Proof. The results follow from Lemma [2] and Theorem [T] and the asymptotic relations given 
in Appendix lAl For the SQB and C1R diffusions: (y>+ , <p~ )(o, x ] < °°i (Vp > fp )(o,x] < 00 
(= oo) if \x S (0,1) (/i > 1) and (<p+, <pf)[ x ,oo) = (^.^)[x,oo) = oo. For the OU diffusion: 
(<Pp ,</ 3 p)(-oo,x] = (^+,Vp)(-oo,a:] = (<£p . </>p ) [*,<») = (^p j Vp ) [x,oo) = oo. Finally we show 
that for the Bessel and confluent families of ^-diffusions (with q 2 > and u e (0,1)) the 
regular boundary F(0+) is killing. By Theorem fl] we need only prove that the regular 
boundary I = of the corresponding process Xf* is killing. Now notice that tpf(x) = 



u p {x) 

properties for the SQB and C1R diffusions from Appendix [X] we have 



? + f f , Res > 0, are fundamental solutions for A( p )-diffusions. Using the asymptotic 



lim tf(*)=0 and lim J_ d JlM = lim ^^^ =ConsM 0. 

Therefore, I = is a killing boundary for x[ p ^ . □ 

Theorem 4. The Bessel, Confluent and OU families of regular diffusions F t — F(X { ), 
with affine drift and monotonic map defined as in Theorem^ conserve the expectation rate, 



i.e., the relation (2.6) holds and hence (e~ bt F t ) t >o is a martingale when the drift parameter 
a — 0, under the following conditions: 

Bessel and Confluent families: c 2 — or c 2 ^ & q 2 ^ & F(0+) = holds; 



OU family: the relation (2.6) holds for every choice of parameters. 



Proof. According to Theorem we need only show that there exists some constant c such 



that the boundary conditions \2.7\ and (2.8) hold when Res > c. Using the respective 



asymptotic relations in Appendix |A| for the SQB, C1R or OU processes, one readily verifies 



that (2.7 1 and (2.8 1 are satisfied for Res > b. □ 
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4.2 Classification of Monotonic Maps 



In addition to the set of monotonic maps F^ 1 ' defined by (1.251, we introduce the following 
four classes of maps: 



Fi 2) (x) 



:(3). 



x) 



Fi 4) (x) 



F^(x) = 



ec 



<pf{x) 
cnp+ +b {x) + c 2 <p~ +b (x) 
<pj(x) 

q%(pp(x) + q2<Pp(x) ' 
cnp+ +b {x) - c 2 ip- +b {x) 
qifp(x) + q2<fip(x) 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



where p, p + b, c, c\, c 2 > 0, e = ±1, and q\, q 2 > with the only exception of F^ 5 ^ for which 
one of <7i, q 2 may be zero. 

As follows from Proposition [2] a function F, given by equation (1.251, (4.51, (4.6|, or 
( 4.7 1 , defines a monotonic map if and only if the Wronskian W{x) — W\u , v 0+ t\(x ) ha s 
the same sign in neighbourhoods of both endpoints / and r. For the map F^ 5 ) in (4.8 1, 
monotonicity follows from Proposition [3] when b > 0. If b < 0, then we again need to 
analyse the asymptotic behaviour of W{x), and hence of W[ip^ ,ip^ +b ](x), as x — > 1+ and 
x — > r— . The asymptotics of such Wronskians are given in Proposition [4] of Appendix [A"| 

The following lemma summarizes the monotonicity properties for all families of maps 
presented so far. 



Lemma 3. Let if be the fundamental soluti ons fo r the S Q B, C IR, or OU diffusion pro- 
cess. The maps F^ 1 ' , I = 1,2, ... ,5, defined in (1.25), (4-5) — (4-8), respectively, are strictly 
monotonic under the following conditions: 

the maps F±' 2 ^ are strictly monotonic for all choices of parameters; 

the map F± is strictly monotonic if b < 0; 

the map F^ is strictly monotonic if b > 0; 

the map is strictly monotonic if and only if b > 0. 

The derivative dFjP /dx, I = 1, 3, 4, has the sign ±e. For cases 2 and 5 we have: sign(dF^' 1 /dx) — 
±esign(6) and sign(dF^/dir) = e. 

Proof. The proof follows directly from Propositions [2] [3] and [4] □ 



5 Conclusions 



By applying the diffusion canonical transformation method which combines a special class 
of variable and measure changes on solvable underlying diffusion processes, this paper has 



G. Campolieti and R. Makarov, Solvable Nonlinear Volatility Diffusion Models 17 



developed new families of exactly solvable multiparameter diffusion models with affine drift 
and nonlinear diffusion coefficient functions. As an application of our general formulation, 
this paper has presented three new main families of time-homogeneous affine diffusions, 
named as Bessel, confluent hypergeometric, and Ornstein-Uhlenbeck families. These pro- 
cesses can be used to model asset prices, interest rates and other financial quantities. In 
particular, analytically exact closed- form expressions for transition densities and state-price 
densities are obtained for these families of nonlinear local volatility models in terms of 
known special functions (i.e., modified Bessel, confluent hypergeometric, hypergeometric). 
The present formulation extends and includes all of the more special families of driftless 
diffusions that were obtained in previous literature with the use of a diffusion canonical 
transformation methodology (e.g., see Subfamilies of these new models encompass 

the CEV model and other related models as special cases, and have been shown to exhibit 
a wide range of implied volatility surfaces with pronounced smiles and skews (see [6j [3 [8] 
for further details). 

By making use of closed-form spectral expansions for the densities of the underlying 
simpler processes, the theory of the present paper can also be used to further derive new 
classes of analytically exact probability densities for first-hitting times and the extrema 
(maximum and minimum) of the new affine drift families of nonlinear local volatility diffusion 
processes. This then also leads to analytical pricing formulas for various lookback and barrier 
options. These applications are the subject of currently related research (see [S]). 
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A Large and Small Argument Asymptotics of the Fun- 
damental Solutions and Wronskian Functions 

Let s, si and s 2 be complex-valued parameters with positive real part. 

A.l The SQB process 

(i) Asymptotic forms of the fundamental solutions: 

2\f2sx/v -2\/2sx/i>a 
+ / \ 1 — f \ 

tpjlx) ~ — == ,„ - -i , -. (p s (x) ~ i J_ ,„ - 1 , . as x — » oo . 
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(ii) Asymptotic forms of the Wronskian functions as x — > 0: 
wr + + u \ /^V s l s 2 V 2(s 2 -si) 

W[<p+,<p-]{x) 



X 



-n-i 



vl ) ^r( M + i)r(/i + 2) ' 

1 (s A V 12 

2 \S2 

II 



2 

x- 2 In ( i) 

r( M )r(i- M ) gf -4 p_i 

4 



if /i = 1 , 

if /i e (0, 1) 



(iii) Asymptotic forms of the Wronskian functions as x — > 00 

1 v^2 - \AI 



4vr (S1S2) 1 / 4 

1 -v/sl + 



2;— /*— 1 e 2( v /sT+ v /3i")v / 25/i/ 



4 (S 1S2 )V4 



4 (S1S2) 1 / 4 



Note that for multiple- valued power function of the form s a with Re s > and a > the 
principle value is used. 

A. 2 The CIR process 

(i) Asymptotic forms of the fundamental solutions: 

~1 ^(x)~^(kx)-» as z^O, 
T(v) 



T(v) 



<p s (x) ~ (kx) ^ as x — > 00 . 



where in Subsections A. 2 and A. 3 we define u = s/Ai. 
(ii) Asymptotic forms of the Wronskian functions as x — > 0: 

(u 2 - 



r( M + i) „ „_! 
r(t*) 

r(m)r(u 3 ) v 1 ^ 

Kr^rtW) x 2m (x) 
K-T^rci-M) / r(v 1 
r(«i)r(t;j) V r ( Ul - 



-2/j 



i> 2 ) 



fi) r(t; 2 -M) 



if /i > 1 , 
if A* = 1 > 

x-v- 1 if/xe(o,i) 



where in Subsections A. 2 and A. 3 we define Vi = Sj/Ai, i = 1,2. 
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(iii) Asymptotic forms of the Wronskian functions as x — > oo: 

w[<p ai ,<p n \w ~ rK)r( W2 ) e {KX) 



— Hi— u 2 — 1 



A. 3 The Ornstein-Uhlenbeck process 

(i) Asymptotic forms of the fundamental solutions: 



^(i)~^( v ^i|r 1 e'« , / J as z^ioo, 
(^J(x) ~ (v / K|a;|) _t ' as x — > =Foo . 



(ii) Asymptotic forms of the Wronskian functions as a; — > — oo: 

rOi)r(u 2 ) 

(iii) Asymptotic forms of the Wronskian functions as x — > oo: 



\V1+V%— 3 gKI 



r(t>i) 



— 1>1 — U2 — 1 



A. 4 Asymptotic Properties of Wronskians of the Fundamental So- 
lutions 

Proposition 4. Let p\ and p 2 be positive real-valued parameters. For the three underlying 
diffusions, namely, the SQB, CIR, and OU diffusions defined in Sections \3.1§3.3\ we have 
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the following limits for the Wronskians of the fundamental solutions: 

W ' ^ J {X) \ sign(p 2 - pi) oo , asx^r-, 

W\<p- <p-](x) - { Sign(pi " P2) °° ' asx ^ l+ > 
^[Vpi'VpjKa:) -00 , as x^l+, 

^[f~pi' L P~t>2\( x ) ~^ ~ 00 ' as 1 ^ r_ 1 / or ^ e CIR and OU 

W[<p+ , ](x) - { > ^ <w * -> r- , /or ifre , 

wii/i constant C > /or i/ie 5"(5B and CIR diffusions and C = /or </ie Of/ diffusion. Here 
sign(x) = +1(— 1) for x > (x < 0) and we define sign(O) = as well as ■ 00 = 0. The 
convergences are monotonic with notation W(x) — > +0 and W(x) — > +00 (or W(x) — > — 
and W{x) — > — 00 j meaning that W{x) > (or W(x) < 0) holds as x approaches an 
endpoint. 

Proof. The limits follow from the asymptotics of the Wronskian functions presented in 
Subsections |A.lf|A.3| where Sj = pi, I — for SQB and CIR, I = —00 for OU and r = 00 for 
all three diffusions. For the CIR diffusion with /1 G (0, 1), the proof for the limiting value of 
W[ifip 1 ,ipp 2 ](x) follows from Proposition |5] □ 

Proposition 5. Let < a < 1 be a real constant. Then: 
(i) S(x;a) — ~ r^jy^ * s an increasing function ofx>0. 
(ii) R(x;a) = r ^ff a ) ^ aw increasing function of x > 0. 

Proof. The first result is proved by 14J. The other is its immediate corollary. Indeed, we 
have 

r(a:) T(x + 1) x-a T(y + a) I y a (y-l) ,y a (y-l) 

-b{y;a)- 



T(x — a) T(x — a + 1) x r(y) y a y — (1 — a) ' y — (1 — a) ' 

where y = 2: + 1 — a>l — a. It is easy to verify that the quotient is an increasing 

function of y. Hence, R{x; a) is a product of two increasing functions. □ 
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Figure 2: Sample local volatility curves for the Conflucnt-W (a) and OU (b) models with 
drift parameters a = 0, b ^ and other parameters calibrated such that <Ti oc (F) = 0.20 at 
F = 100. 



